Structural properties and optical response of Na clusters in Ne, Ar, and Kr matrices 
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We discuss the structural properties and optical response of a small Na cluster inside rare gas 
(RG) matrices of Ne, Ar, or Kr atoms. The mixed systems are described with a hierarchical model, 
treating the cluster at a quantum mechanical level and the matrix atoms classically in terms of 
their positions and polarizations. We pay special attention to the differences caused by the different 
matrix types. These differences can be explained by the interplay of core repulsion and dipole 
attraction in the interaction between the cluster electrons and the RG atoms. 

PACS numbers: 36.40.Gk,36.40. Jn,36.40.Qv,36.40.Sx 
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I. INTRODUCTION 

The study of the properties of clusters embedded in a 
matrix or deposited on a substrate has motivated many 
researches since several years 1 ^. This setup becomes in- 
creasingly important because the embedding/depositing 
simplifies the experimental handling and because it is a 
generic test case for composite materials. As in free clus- 
ters, the doorway to (laser induced) cluster dynamics is 
the optical response, especially in metal clusters where 
the Mie plasmon dominates the optical properties 3 -. Let 
us mention as examples the systematics of optical re- 
sponse in large noble-metal clusters^iSil and its depen- 
dence on the environment^. The study of the optical 
response thus constitutes the key for understanding the 
response of clusters to electromagnetic probes and it also 
serves as a powerful tool for analyzing the underlying 
cluster structure. 

The case of inert environments is especially interesting 
because it implies only moderate perturbations of cluster 
properties. One can thus benefit from the well defined 
conditions from the surrounding system and still access 
predominantly the cluster properties. But the theoretical 
modeling becomes much more involved and the develop- 
ment of reliable as well as inexpensive approaches is still 
a timely task, in particular what truly dynamical appli- 
cations arc concerned 9 - ' 10 ' 11 . Nonetheless, fully detailed 
calculations have been undertaken where details count, 
e.g. for the structure of small Na clusters on NaCli^ 
or the deposit dynamics of Pd clusters on a MgO sub- 
strate^. But the expense for a fully fledged quantum 
simulation grows huge. These subtle models are hardly 
extendable to truly dynamical situations, to larger clus- 
ters or substrates, and to systematic explorations for 
broad variations of conditions. This holds true not only 



for clusters on substrates, but for all composite systems. 
Thus there exists a great manifold of approximations 
which aim at an affordable compromise between reliabil- 
ity and expense. One route, for example, keeps all con- 
stituents at the same level, but simplifies the description 
in terms of a microscopically founded tight-binding ap- 
proac h 14 i 15 i 16 i 17 i 18 i 19 . At the other extreme, one can con- 
sider all degrees of freedom as classical and perform pure 
molecular dynamics, as e.g. the deposition dynamics of 
Cu clusters on metal^ - or Ar— surface, and of Al or Au 
clusters on SiO^. In between, one can take advantage 
of the very different importance or activity within the 
composite and thus develop a hierarchical modeling us- 
ing various levels of approximation for the different sub- 
systems. Such approaches are widely used in quantum 
chemistry, often called quantum-mechanical-molecular- 
mechanical (QM/MM) model. They have been applied 
for instance to chromophores in bio-molecule a 23 ' 24 , sur- 
face physic s 25 ' 26 , materials physic s 27 ' 28 ' 29 ' 30 , embedded 
molecules 3 ^ and ion channels of cell membranes^ 2 .. A 
variant of mixed modeling was applied in the case of a 
Cs atom in He environment, the latter also considered as 
a quantum system^ 3 -. 

We are dealing here with Na clusters in rare-gas en- 
vironments (Ne, Ar, Kr). The large difference between 
cluster metals (reactive) and rare gas (inert) naturally 
suggests a hierarchical model where the substrate atoms 
are handled at a lower level of description, as classical 
particles but with a dynamical polarizability. Taking up 
previous developments fro m 34 ' 35 , we have developed, in 
the spirit of QM/MM approaches, a hierarchical model 
for Na clusters in contact with Ar— and applied it to 
structure, optical response^ 7 - and to non-linear dynam- 
ics of embedded clusters^ 8 -. It is the aim of this paper 
to present a generalization to other types of rare gases 
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(RG), namely Kr and Ne, and a first comparative study 
of the effects of different environments on structure and 
optical response. As test cases, we consider a Nag cluster 
embedded in RG clusters of various sizes. Strictly speak- 
ing, they are mixed clusters and properties which depend 
on the size of the RG system are specific to mixed clus- 
ters. But we use the mixed systems mainly as model for 
Na clusters embedded in a matrix and we thus consider 
rather large systems. Henceforth, we will use the notion 
"matrix" for the RG surroundings. 



II. MODEL 

The model has been introduced and presented in detail 
ir*22. However for sake of completeness, we recall in this 
section the ingredients and a few relevant formulae. 

The degrees of freedom of the model are the wave- 
functions of valence electrons of the metal cluster, 
{(p n (r),n = l...iV c i}, the coordinates of the cluster's 
Na+ ion cores, {R/,/ = l...iV ion }, of the Ar atoms 
cores Ar Q+ ), {R a ,a = I...AA r }, and of the Ar va- 
lence clouds, {R' a ,a = f...AAr}- From the given total 
energy, the corresponding equations of motion are de- 
rived in a standard manner by variation. This leads to 
the (time-dependent) Kohn-Sham equations for the one 
single-particle wavefunctions (p n (r) of the cluster elec- 
trons, and Hamiltonian equations of motion for the other 
three degrees of freedom, thus treated by classical mole- 
cular dynamics (MD). For the valence cluster electrons, 
we use a density functional theory at the level of the 
time-dependent local-density approximation (TDLDA), 
augmented with an average-density sclf-intcraction cor- 
rection (ADSIC)^. The density of these electrons is 
given naturally as defined in mean-field theories and 
reads p e i(r) = J2n l</?ri( r )| 2 - A RG atom is described by 

I 



two constituents with opposite charge, positive RG core 
and negative RG valence cloud, which allows a correct 
description of polarization dynamics. In order to avoid 
singularities, we associate a smooth (Gaussian) charge 
distribution to both constituents having width cr RG of 
the order of the p shell "size" in RG atoms, in the spirit 
o£i: 
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The corresponding Coulomb potential exerted by the RG 
atoms is related to the charge distribution (fT|) by the 
Poisson equation, and reads: 



erf(|r-R a |/cr RG ) 



|r-RJ 
erf(|r-R^|/a RG ) 
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(2) 



where erf(r) = dx e x stands for the error func- 

tion. As for the Na + ions, their dynamical polarizability 
is neglected and we treat them simply as charged point 
particles. 

The total energy of the system is composed as: 
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The energy of the Na cluster -ENaciuster consists out of 
TDLDA (with SIC) for the electrons, MD for ions, and 
a coupling of both by soft, local pseudo-potentials, for 
details se o n i 41 ' 42 . The RG system and its coupling to 
the cluster are described by 
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, (4) 
(5) 

(6) 



The Van der Waals interaction between cluster elec- 
trons and RG dipoles is written in Eq. ([6]) as a correlation 
from the dipole excitation in the RG atom coupled with 
a dipole excitation in the cluster, using the regularized 
dipole operator f a corresponding to the smoothencd RG 
charge distributions^^. 

The interaction of one RG atom with the other con- 



stituents (RG atoms, Na + ions, cluster electrons) results 
from the balance between a strong repulsive core poten- 
tial that falls off exponentially and an equally strong at- 
traction from dipole polarizability. The (most important) 
polarization potential is described by a valence electron 
cloud oscillating against the RG core ion. Its parameters 
are the effective charge of valence cloud Q, the effective 
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mass of valence cloud mAr = Qrn^ the restoring force 
for dipoles /org, an d the width of the core and valence 
clouds org- The Q and &rg ar c adjusted to reproduce 
experimental data on dynamical polarizability od(cj) of 
the RG atom at low frequencies, namely the static limit 
au{oJ = 0) and the second derivative of a'^(J' = 0). 
The width org is determined consistently such that the 
restoring force from the folded Coulomb force (for small 
displacements) reproduces the spring constant /crg- 

The short range repulsion is provided by the vari- 
ous core potentials. For the RG-RG core interaction in 
Eq. ([¥]), we employ a Lennard- Jones type potential with 
parameters reproducing binding properties of bulk RG : 



Vr 



(core) 



RG.RG v 



(R) = e 2 A RG {Rrg/R) 12 - (Rrg/R) 



(7) 



The Na-RG core potential V^ a rg in Eq. ([5]) is chosen 
according tc4£, within properly avoiding double counting 
of the dipole polarization-potential, hence the following 
form : 
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The parameters are taken from literature for Na-Nei 4 -, 
Na-AriM£ an d Na-Kri&. The pseudo-potential W e i,RG 
in Eq. j5|) for the electron-RG core repulsion has been 
modeled according to the proposal oS^ : 
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with a final slight adjustment to the properties of a Na- 
RG molecule (bond length, binding energy, and optical 
excitation spectrum). Values of atomic and dimcr prop- 
erties used are reported in table |TJ 



RG 


RG-Atom 


Na-RG 


RG bulk 




org [a 3 ] 


IP [Ry] 


do [ao] 


E () [mRy] 


r s [ao] 


E coh [mRy] 


Ne 


2.67 


1.585 


10.01 


0.0746 


5.915 


-0.0272 


Ar 


11.08 


1.158 


9.47 


0.3793 


7.086 


-0.1088 


Kr 


16.79 


1.029 


9.29 


0.6238 


7.540 


-0.1497 



TABLE I: Properties of the RG atoms and Na-RG dimers 
for Ne, Ar, and Kr which were used for the fine-tuning of the 
model. 



III. MODUS OPERANDI 

The numerical solution proceeds with standard meth- 
ods as described in detail ini2. The TDLDA equations 



for the cluster electrons are solved on a grid in coordi- 
nate space, using a time-splitting method for the propa- 
gation and accelerated gradient iterations for the station- 
ary solution. We furthermore employ the cylindrically- 
averaged pseudo-potential scheme (CAPS) as an approx- 
imation for the electron o 47 ' 48 , which is justified for the 
chosen embedded Na 8 cluster. We have checked that a 
2D calculation with CAPS and a full 3D treatment of the 
valence electron wavefunctions both give almost identical 
optical responses. The Na + ions as well as the RG atoms 
are treated in full 3D. The dynamics of the Na electrons 
is coupled to the response of the RG dipoles. However, 
the ionic and atomic positions can safely be frozen for 
the present study as we focus on exploring the optical 
response of the embedded clusters. 

To find an optimal Na+RG configuration, one starts 
with a fee RG crystal, cuts from that a given number of 
closed shells, and cools the resulting configuration for a 
pure RG cluster. One then carves a cavity of 13 atoms 
(Ar, Kr) or 19 atoms (Ne) from the center and places 
the Nas cluster into it. This mixed configuration is re- 
optimized by means of successively cooled molecular dy- 
namics for the ions and atoms coupled to the stationary 
solution for the cluster electrons. 

The stationary solution of the equations of motion pro- 
vides the ground state of the mixed system and consti- 
tutes the initial condition for further dynamical calcu- 
lations. One can then compute several observables to 
analyze both the statics (structural properties) and the 
dynamics. A global measure for ionic and electronic clus- 
ter structure arc the r.m.s. radii. 
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where (. . .)/ = J^i ■ ■ ■ ari<1 (■ • -}e = J d 3 r p e (r) . . . Note 
that these quantities may also be used for characteriz- 
ing dynamics, although they take interest mostly on very 
long times when ions and atoms actually move. This lat- 
ter aspect is not directly addressed here and we shall thus 
use them only as static quantities. We also evaluate the 
insertion energy Ei ns which is defined as 



-Ems — -Etot(NasRG N ) 
— -E(RGn+p) 



f E(RG P ) 
£(Na 8 ) 



(11) 



with p = 13 in the case of Ar and Kr, and p = 19 for 
Ne. A further useful energetic observable is the ionization 
potential (IP), that is the energy required to remove one 
electron from the metal cluster. We compute it as the 
single-particle energy of the least bound electron which 
is a reliable measure in ADSIG^. 

At the side of truly dynamical properties, as already 
emphasized, one should remind the especially important 
role played by the optical response. This observable is 
computed in an explicitely dynamical way. The dynamics 
is initiated by an instantaneous dipole boost of the clus- 
ter electrons. The optical response is then obtained by 
spectral analysis of the emerging time-dependent dipole 
signal following the strategy proposed i n 49 i 50 i 51 . 
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IV. RESULTS AND DISCUSSION 

Figure [T] shows the RG structure in terms of radial dis- 
tributions of atoms. The effect of embedding is visualized 
by comparing the pure RG system (dotted) with the RG 
distribution around the Nas cluster (full lines). The dis- 
tributions line up nicely in radial shells. For pure RG 
clusters, they remain very close to the radial shells of the 
bulk fee structure (not shown here, for the case of Ar, 
see^ 7 -). The overall scale is basically given by the bulk 
Wigner-Seitz radius r s given in table [J Note that Ne 
has a much smaller r s resulting in denser packing as seen 
in figure [TJ Carving of the cavity and insertion of Nag 
has only small effect in the Ar and Kr environment and 
mainly for a few inner shells. For Ne, however, we see a 
stronger perturbation which spreads over all atoms. The 
examples in figure [1] concern rather larger RG systems 
probably close to bulk. The embedding effects increase 
with decreasing number of RG shells, remaining small 
throughout for Kr and Ar, but soon destroying any clear 
shell structure for Ne^. The reason is that Ne is much 
less bound than Ar or Kr as can be read of from the 
cohesion energy in table [H and already these seemingly 
more robust materials are weakly bound. 

The impact of the RG environment on metal cluster 
properties is analyzed in figure [H which shows global ob- 
servables of the ground state configurations of Nas m 
the various matrices, as a function of matrix size. The 
insertion energies in the uppermost panel indicate that 
Ne cannot finally capture Nas inside while Ar, and even 
more so Kr, provide robust environments for embedded 
metal clusters. The binding increases slowly with matrix 
size, except for Ne where the polarizability is too weak to 
accumulate sufficient long range attraction. The IP (sec- 
ond row in figure[2]) makes a jump down from free Nas to 
the embedded cases and then stabilizes with a few fluc- 
tuations for small matrices and a faint further decrease 
for larger ones. The energies are at the same scale for 
all RG types. The sudden drop from free to embedded 
for Ar and Kr is due to the core repulsion from the first 
RG shell exerted on the cluster electrons. Adding further 
shells acts only indirectly by compressing the whole ma- 
trix and thus bringing the innermost shell slightly closer 
to the cluster. Radii in Nas a re shown in the two lower 
rows. They vary very little in general. There remain 
interesting differences in detail. The trends with matrix 
size are the same for electrons and ions. However the 
step from free to embedded is much different to the ex- 
tent that electrons are more compressed when embedded 
which is a visible effect from core repulsion, similar as the 
jump in IP. The radii decrease slightly with matrix size 
for Ne and increase for Ar or Kr. This indicates that core 
repulsion prevails in Ne while dipolc attraction becomes 
more effective in Ar and Kr. 

We finally present in figure[3]the trends of the Nas plas- 
mon resonance peak. The changes are generally small, 
at an absolute scale at limits of our modeling (which 
we estimate to about 0.1 cV uncertain). The relative 



trends, however, can be taken at smaller energy scale 
and these carry several interesting aspects. As previ- 
ously discussed^!, the position of the peak results from 
a subtle cancellation between core repulsion and polar- 
ization effects. The step from free to embedded clusters 
first produces a blue shift because the cluster electrons 
feel the core repulsion from the first layer of RG atoms. 
The plasmon peak moves slowly back to red with increas- 
ing system size because each new RG shell adds to the 
long-range and attractive polarization potential. Polar- 
izability is also the key to the trend with RG material. 
It increases with atomic number (see column 2 of table 
|T|. The cancellation mentioned above is thus more effec- 
tive in Kr than in Ar, and in Ar more than in Ne. This 
explains the steady decrease of the peak position from 
Ne to Kr and the growth of slope with atomic number. 
In fact, Ne shows no significant slope to red at all. This 
happens because added atoms compress slightly the in- 
nermost shell (see shrinking radii with increasing system 
size in figure [5]) and enhance core repulsion which, in 
turn, compensates the growth of polarization effects. Ar 
and Kr experience no such compression and have anyway 
the stronger polarizability. 

Experimental data for Na clusters embedded in rare 
gas material are not yet available. A direct compari- 
son has thus to be postponed. But there exist already 
some data on the optical response in the somehow sim- 
ilar combination of Ag clusters in RG material: the 
Lausanne group has studied, among others, Ags@Ar— , 
Ag 7 @Ar,Kr,Xe^, and Ag 7 @Ar,Ne^ in large RG matri- 
ces and, the Rostock group Ag 8 covered by small layers 
of Ne, Ar, Kr and Xe (of sizes between 4 and 135 RG 
atoms) all immersed in a He droplet^. 

Table |TT] summarizes the results in terms of shifts of 
the plasmon peak position with changing RG material. 
A direct comparison with free neutral Ag clusters is not 
available because these are hard to handle experimen- 
tally. As for experiments performed in£, Ags is embed- 
ded in very small RG matrices, and the whole system 
itself embedded in a He droplet, for reasons of better 
handling. It is usually claimed that the helium envi- 
ronment interacts very faintly with the embedded sys- 
tem and can then be considered practically as a vacuum. 
However experimental dat a 56 ' 57 and DFT calculations^ 
show that the presence of He around Cs atoms produces 
a blue shift of the plasmon peak. It is thus likely that 
some small blue shift exists also for Ag clusters directly 
in He droplets. Therefore a comparison with truly free 
clusters is excluded. On the contrary, the influence of 
He around the RG layers in the experiments o£^ is most 
probably negligible. The effect of the He droplet is in- 
deed strongly shielded by the RG layers both spatially 
and energetically. Remind that a typical RG-He bond 
has a length of about 6-7 and an energy of a few 
meV— , comparable to the metal-RG bonds^^^ (see 
fourth column of table fl} . But once coated with the RG 
layers, the metal cluster lies typically twice farther away 
from the He droplet than in the case without RG layers, 
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FIG. 1: Radial distribution of RG atoms in Na8Ne428, NagAr434, and NasKr434 (full lines) compared with the distributions for 
the pure RG cluster with 447 atoms before carving the cavity (dashed lines). 
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FIG. 2: Ground state observables (energies, radii) for Na8 embedded in Ne, Ar, and Kr matrices of different size. 
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He drop -> Nei 35 : 0.009 


Nc -» Ar: -0.1 






He drop -> Ar: -0.03 


He drop -> Ar 24 : -0.01 


Ar -> Kr: -0.06 / - 0.09 








He drop Kr 25 : -0.04 
He drop -> Xei 5 : -0.08 



TABLE II: Shift of plasmon peak position, in eV, going from one environment to another, as indicated. A plus (minus) sign 
stands for a blue (red) shift. 
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FIG. 3: Plasmon resonance energies for Nas in matrices of 
different size and RG material. 



by the ratio of Wigncr-Seitz radii, that is 3/4, and by 
reoptimizing the RG positions before the calculations of 
the optical response of this pseudo Ag 8 . We find the 
same trends as with Nag in figure [31 but indeed reduced 
by a factor of three. 

Besides, previous theoretical calculations were per- 
formed within TDLDA and a jcllium model, coupling to 
the RG materials in terms of a static dielectric medium 
with dielectric constant e for different metal clusters, 
namely K-2I, Na and Al2£, and Ag22. They always yield 
red shifts which, moreover, increase with increasing e. 
Indeed, the ingredients added in these models can only 
produce long-range polarization effects and thus can only 
give a red shift. Our model contains as new component, 
the RG short-range and repulsive core potentials which 
generate a blue shift of the plasmon peak. A clear exper- 
imental assessment would yet require a comparison with 
truly free neutral metal clusters. 



whence a vanishingly small residual interaction between 
the metal cluster and the He droplet. Thus the rela- 
tive shifts that we observe fby changing RG material are 
most probably reliable. Finally, a word of caution is in 
order. The above mentioned experimental measurements 
and our calculations on embedded Na clusters^ report 
a broadened, often even fragmented, peak with a width 
somewhat larger than the shifts we are looking at. The 
comparison is thus at the edge of experimental and the- 
oretical resolution. 

The data in table [UJ agree with our theoretical results 
in that all shifts are very small. Inert environment turns 
out to be indeed inert with respect to the plasmon peak 
position. Looking in more detail at the relative shifts, 
one can read for the step from Ne to Ar a red shift of 
0.1 eV for Ag7 in large matrices^ and of 0.02 eV for 
Agg in small matrices^. The step from Ar to Kr yields 
0.06-0.09 eV for Ag7 in large systems^ and 0.03 eV for 
Agg in small systems^. The experiments with small RG 
layers in a He droplet thus yield generally smaller red 
shifts. This also holds for the step from pure He to Ar 
environment. This is probably explained by the size of 
the matrix, as we see from our results in figure [3] a slow 
but steady move towards red (however, systematic er- 
rors by comparing two very different experiments cannot 
safely be excluded). Our results for embedded Na clus- 
ters shown in figure [3] show a red shift of 0.1-0.2 eV for 
the step from Ne to Ar and of 0.15-0.2 eV for Ar to Kr, 
both growing with increasing system size. They confirm 
all trends seen in experiment but are generally half an 
order of magnitude larger. This is probably due to the 
smaller Wigner-Seitz radius of Ag (3 ao instead of 4 ao 
for Na) which means that Ag structures are much more 
compact than Na ones and thus couple less strongly to 
the matrix because both fill the same RG cavity. And 
this, in turn, produces smaller shifts. In order to check 
that argument, we have simulated an embedded "Agg" 
cluster simply by rescaling the ionic positions of the Nag 



V. CONCLUSIONS 

We have discussed in this paper static properties and 
optical response of a small Na cluster embedded in var- 
ious rare gas (RG) matrices. For this purpose, we used 
a recently introduced hierarchical approach combining 
a fully detailed quantum-mechanical description of the 
cluster with a classical modeling for the RG environment 
and its interactions with the cluster. We have studied 
effects from embedding in RG environment with up to 
434 atoms with particular emphasis on the change with 
RG material. Various observables were considered. The 
insertion energy yields stable embedding for Ar and Kr 
but not for Ne within the considered system sizes. The 
IP behaves the same in all three materials: It drops from 
free to embedded and stays nearly constant for all ma- 
trix sizes. The electronic and ionic radii of the Na cluster 
change very little; the most noteworthy effect is here a 
slight compression of the electron cloud through embed- 
ding. What optical response is concerned, we studied the 
effect of the RG matrices on the position of the surface 
plasmon peak in the metal cluster. The net shift from 
free to embedded clusters is very small due to a near 
cancellation of the blue shift from core repulsion with 
the red shift from dipole polarization in the interaction 
with the RG atoms. The polarization increases with the 
RG atomic weight and thus final peak position goes from 
blue shift to red shift on the way from Ne over Ar to 
Kr. The polarization effect also increases with the size of 
the RG matrix which produces a slow and steady trend 
to red with increasing system size. Comparison with ex- 
perimental data on embedded Ag clusters confirms these 
trends and orders of magnitude. Having explored the 
basic observables of structure and optical response, the 
way is open to applications in truly dynamical scenarios. 
The large difference in atomic masses of the RG will then 
play a crucial role and yield interesting effects. Work in 
that direction is under progress. 
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